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ABSTRACT 

We begin with a review of a well known class of networks, Classical Bayesian (CB) 
nets (also called causal probabilistic nets by some). Given a situation which includes 
randomness, CB nets are used to calculate the probabilities of various hypotheses 
about the situation, conditioned on the available evidence. We introduce a new class 
of networks, which we call Quantum Bayesian (QB) nets, that generalize CB nets to 
the quantum mechanical regime. We explain how to use QB nets to calculate quantum 
mechanical conditional probabilities (in case of either sharp or fuzzy observations) , 
and discuss the connection of QB nets to Feynman Path integrals. We give examples of 
QB nets that involve a single spin-| particle passing through a configuration of two or 
three Stern-Gerlach magnets. For the examples given, we present the numerical vàlues 
of various conditional probabilities, as calculated by a general computer program 
especially written for this purpose. 
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1. INTRODUCTION 

The artificial intelligence and expert systems literature contains a large number of 
articles discussing the theory and application of Classical Bayesian (CB) nets. 

At least one software package|| that implements this theory is commercially available. 
Amazingly, the physics community, to whom this paper is mainly addressed, seems 
not to have discovered or used CB nets yet. Therefore, we begin this paper with a 
review of CB nets. However, the real purpose of this paper is to introduce a new 
class of nets, which we shall call Quantum Bayesian (QB) nets, that generalize CB 
nets to the quantum mechanical regime. In this paper, to illustrate QB nets, we use 
them to predict conditional probabilities for experiments comprising combinations 
of two or three Stern-Gerlach magnets. In future papers, we will use QB nets to 
analyze quantum optical experiments and other physical situations. Earlier workers 
have devised networks for quantum mechanics ( for eg., Feynman diagrams, and the 
"trajectory graphs" of Ref.|TtJ), but their nets differ substantially from ours. 



Henceforth, we will underline random variables. We will write P(x = x) for the 
probability that the random variable x assumes the particular value xfnj. Sometimes, 
if there is no danger of confusion, we will write P(x) rather than P(x = x). Similarly, 
we will often write P(y\x) instead of P(y — y\x — x). P(y — y\x — x), the conditional 
probability that y assumes the value y given that x assumes the value x, is defined 
by 

P(y = y, x = x) 

P{y = y\x = x) = r . (1.1) 

- P(x = X) 

Suppose h is a random variable that can take on as vàlues each of a number of 
hypotheses, and e is a random variable representing the observations that constitute 
the available evidence. We often want to calculate the posterior probability P(h\e) in 
terms of the prior probability P{h). To do this, one may use Bayes's rule, which says 

P(h\e)= (12) 

Equation(1.2) is a simple consequence of Eq.(l.l). The subject of CB nets may be 
viewed as an extension of Bayes's rule. 

A CB net has two parts: a diagram consisting of nodes with arrows connecting 
some pairs of these nodes, and a collection of probabilities, one per node. For example, 
the digital circuits of electrical engineering can be modelled as CB nets if each NAND 
(or, alternatively, each AND, OR and NOT gate) is replaced by a node, and each 
connecting cable is replaced by an arrow pointing in the direction of current flow. 
Of course, in the usual digital circuits, the NAND gates are deterministic, their 
output being a deterministic function of their inputs. More generally, if some of 
the NAND gates were to act erratically, and if the input sources to the circuit also 
acted in a random fashion, then such a circuit could also be modelled by a CB 
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net. In this non-deterministic CB net, the signal flowing in each wire would have 
a certain probability of being and of being 1. If the signal in a certain wire were 
measured so that one knew that it was definitely 1, then one would want to revise 
the probability distributions for the signals in all other wires so as to reflect the 
new evidence. Hence, conditional probabilities and Bayes's rule arise naturally when 
considering probabilistic nets. 

As we shall see, a major difference between CB nets and QB nets is that for 
QB nets, one assigns complex amplitudes rather than probabilities to each node. 

2. THEORY OF CB NETS 

In this section, we will review the simple theory of CB nets. The next section will 
present some examples. 

We call a graph (or a diagram or an architecture) a collection of nodes with 
arrows connecting some pairs of these nodes. The arrows of the graph must satisfy 
certain constraints that will be specified below. We call a labelled graph a graph whose 
nodes are labelled. A CB net consists of two parts: a labelled graph with each node 
labelled by a random variable, and a collection of node matrices, one matrix for each 
node. These two parts must satisfy certain constraints that will be specified below. 

We define two kinds of arrows: infernal arrows are those that have a starting 
node and a different ending one; external arrows are those that have a starting node 
but no ending one. We define two types of nodes: external nodes are those that have 
a single external arrow leaving it, and internal nodes are those that have one or more 
infernal arrows leaving it. It is also common to use the terms root node or prior 
probability node for a node which has no incoming arrows, only outgoing ones. We 
restrict all nodes of a graph to be either internal or external. Hence, no nodes may 
have both an external and one or more internal arrows leaving it. 

We define each node of a CB net to represent a numerical value, and the 
whole net to represent the product of all these node vàlues. We assign a numer- 
ical value to each node as follows. First, we assign a random variable to each 
node. Suppose the random variables assigned to the N nodes are x_ ^ , , * * * , x m . 
Define Z N = {1, 2, • • • , N}. For any finite set S, let \S\ be the number of ele- 
ments in S. If S = {£4, k 2 , ■ ■ ■ , k\ S \} C Z N , define (x.) s = (x kl , x k2 , ■ ■ ■ , x k]s] ) and 
(x-)s — x k2 , ' ' ' ) Sometimes, we also abbreviate (x.)z N (he., the vector 

that includes all the possible Xj components) by just x., and (x.)z N by just x. . For 
j G Zn, we imagine node x_j to lie in state Xj (See Fig.l). We also imagine all arrows 
leaving the node Xj to lie in state Xj, and thus we label all of them Xj. At this point 
we've shown how to label each arrow in the graph by x k for some k G Z^. Define Sj 
to be the set of all k such that an arrow labelled x k enters node Xj. Now we assign 
a value P[xj\(x.) Sj \ to node Xj. P[xj|(x.)s.] is what we referred to earlier as a node 
matrix; Xj is the matrix's row index and (x.)s j is its column index. As the notation 
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suggests, we assume that the vàlues P[xj\(x.) Sj \ are conditional probabilities; i.e., 
that they satisfy 

P[^|(x.)5j>0, (2.1) 

E^K^-hi, (2.2) 

Xj 

where the sum in Eq.(2.2) is over all states Xj that the random variable x_j can assume, 
and where Eqs. (2.1) and (2.2) must be satisfied for all j G Zj^ and for all possible 
vàlues of the vector (x.)s- of random variables. The left hand side of Eq.(2.2) is just 
the sum over the entries along a column of a node matrix. The CB net is taken to 
represent the product of all the probabilities P[xj\(x.) Sj ] for j E Z N . This product is 
a function P(x.) of the current states xi, x 2 , • ■ • , x^ of the nodes. Thus, 

P(x.)= n P[xj\(x.) Sj ] . (2.3) 

We require that 

E^.) = l, (2-4) 

X. 

as expected for a joint probability distribution of the random variables x_y, x 2 , • • • , x N . 

Next, to illustrate the CB net concepts just presented, we will discuss all 
possible CB nets with two or three nodes. 

For 2 nodes labelled by random variables x and y, there are only 2 possible 
labelled graphs, depending on whether the arrow points from xtoy (Fig.2a) or vice 
versa (Fig.2b). 

Figure 2a is diagrammatic notation for the following equation: 

P(x,y) = P(y\x)P(x) . (2.5) 

Notice that each probability factor on the right hand side of Eq.(2.5) is represented by 
a node in Fig.2a. Notice also that the prior probability P(x), since it has no conditions 
placed on it, is portrayed in Fig.2a by a node which has no arrows pointing into it. 
The probability P(y\x), on the other hand, depends on the value x of node x, and 
thus it is represented by a node with an arrow labelled x pointing into it. According 
to Eq.(l.l), Eq.(2.5) is a tautology (i.e., a statement that is always true). Thus, the 
net of Fig.2a may represent any probability distribution of two random variables x 
and y. 

Figure 2b is diagrammatic notation for the following equation: 

P(x,y) = P(x\y)P(y) . (2.6) 
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Again, this equation is a tautology. Thus, like the net of Figs.2a, the net of Fig.2b 
may represent an arbitrary probability distribution of two random variables x and y. 
(A collorary is that two CB nets with different labelled graphs may still represent the 
same probability distribution. ) 

Figure 3a is diagrammatic notation for 

P(x,y,z) = P(z\x,y)P(y)P(x). (2.7) 

Contrary to Eqs.(2.5) and (2.6), Eq.(2.7) does not represent a tautology: not all 
probability distributions of three random variables x, y and z must satisfy Eq.(2.7). 
In fact, summing both sides of the last equation over the vàlues of z yields 

P(x,y) = P(y)P(x) , (2.8) 

i.e., x and y are independent. Thus, the net of Fig.3a represents two independent 
random variables x and y. 

Figure 3b is diagrammatic notation for the following equation: 

P(x, y, z) = P(z\x)P(y\x)P{x) . (2.9) 
Dividing both sides of Eq.(2.9) by P(x) yields: 

P(y, z\x) = P(z\x)P(y\x) . (2.10) 

Even though y and z are not necessarily independent, they are conditionally inde- 
pendent, i.e., they are independent at fixed x — x. Thus, the net of Fig.3b represents 
two conditionally independent random variables y and z. 

Figure 3c is diagrammatic notation for the equation: 

P(x,y,z) = P(z\y)P(y\x)P(x). (2.11) 

Random variables x_i,x 2 , ■ ■ ■ ,x N form an N step Markov chain (N may be infinite) 
if P(x n+ i\x n ,x n - 1 , ■ ■ ■ ,zi) = P(x n+ i\x n ) for n = 1, 2, • • • , N - 1, i.e., if the (n+l)th 
step x n+ \ depends only on the step x n immediately preceding it. Clearly, the random 
variables x, y and z of Fig.3c represent a 3 step Markov chain. 
Figure 3d is diagrammatic notation for the equation: 

P(x, y, z) = P(z\y, x)P(y\x)P(x) . (2.12) 

Replacing the conditional probabilities on the right hand side of Eq.(2.12) by their 
definitions in terms of probabilities without conditions, we obtain P(x, y, z) for the 
right hand side. Thus, Eq.(2.12) is a tautology, and the net of Fig.3d can represent 
any probability distribution of three random variables x, y and z. 

The graph of Fig.3d is acyclic. That is, it does not contain any cycles (a cycle 
is a closed path of arrows with the arrows all pointing in the same sense). On the 
other hand, the graph of Fig.4 is paracyclic (i.e., it contains at least one cycle). In 
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fact, the whole graph of Fig.4 is a cycle. The net of Fig.4 shall be forbidden because, 
if one sums over its free indices (x, y, z), one does not always obtain unity, as must be 
the case if the net is to represent a probability distribution P(x,y,z). For example, 
if we assume that x, y and z can assume only states and 1, and if P{y\x) = 5 VjX , 
P(x\z) = 5 XtZ and P(z\y) = 5 Zt y, where S x>y is the Kronecker delta function, then Fig.4 
adds up to J2 x ,y,z 5x,z5z,y5y, x — 2. The acyclic net of Fig.3d does add up to one, as 
proven in diagrammatic notation in Fig.5. In this figure, summation over the states 
of an arrow is indicated by giving the arrow a double shaft. Note that in Fig.5, we 
first add over the index z of the external arrow. This produces a new external arrow, 
and we add over its index y, and so on. Each time we add over the index of the 
current external arrow until finally we get unity. One cannot follow this summation 
procedure to show that the net of Fig.4 adds up to one. Indeed, Fig.4 has no external 
arrow to begin the procedure with. 

Certain aspects of the preceding discussion of nets with two or three nodes 
can be generalized to any number iV of nodes. 

For any number of nodes, we will say that a graph is fully connected if it 
is acyclic, and if everyone of its nodes is connected to all other nodes by either 
incoming or outgoing arrows. We will say that a net is fully connected if its graph 
is fully connected. Any fully connected CB net represents a completely general joint 
probability distribution of the random variables labelling its nodes. Indeed, it is 
always possible to label the nodes of an N node fully connected CB net so that the 
net represents the right hand side of the following equation: 

P(X.) = P(X N \X N _ 1 ,X N _ 2 , • • • , Xi)P{x N _i\x N _ 2 , %N-3 • • • , Xi) • • • P(x 2 \x 1 )P(x 1 ) . 

(2-13) 

And this last equation is a tautology. To label the nodes of a fully connected net so 
that Eq.(2.13) applies to it, one proceeds as follows. There always exists exactly one 
external node. This is why The external node, call it xn, must exist. Otherwise, all 
nodes would have at least one outgoing infernal arrow. Then one could start from 
any node and travel along one of its outgoing infernal arrows to reach another node, 
and so on, until one came back to a previously visited node. Therefore, the graph 
would not be acyclic. The external node xn is unique because, since all other nodes 
have outgoing arrows that point to xn, all other nodes must be infernal. Remove 
Xn from the graph. For the same reasons as before, the resulting diminished graph 
contains a unique external node. Call the latter node xn-i- Continue removing nodes 
in this way until all the nodes are labelled xi, x 2 , ■ • ■ , x^. Since Eq.(2.13) suggests 
that Xj occurs after Xj-i for j = 2, 3, • • • , N, we call this node labelling (ordering) 
the chronological labelling of the graph. Two colloraries of the preceding proof are 
that fully connected graphs have a single external node, and that all fully connected 
iV-node labelled graphs are identical once they are relabelled in the fashion described 
above. Figure 6a shows the fully connected four node graph with its chronological 
labelling. By deforming Fig.6a into a topologically equivalent diagram, one obtains 
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Fig.6b, a more "stylized" version of the same thing. Fig.6b might make more clear 
to some readers how the arrows of a fully-connected graph are organized. 

Note that one can relabel any graph chronologically, even if it isn't fully con- 
nected. Indeed, given any graph G, one may add arrows to it to form a fully connected 
graph G. Call G a completion of G. Sometimes it is possible to add arrows to G in 
two different ways and arrive at two different completions. Hence completions are not 
always unique. Any completion of a graph G may be labelled chronologically follow- 
ing the procedure described above. This, in turn, gives the graph G a chronological 
ordering, albeit, not a necessarily unique one. Suppose that x and y are two nodes 
in a graph G. We say that x precedes y and write x < y if , for any completion of G, 
there exist integers j and f with j < f so that x = Xj and y = Xj,. We say that x and 
y are concurrent and write x ~ y if x precedes y in some completions but y precedes 
x in others. And we say that x succeeds or follows y and write x > y if y < x. 

Call a CB pre-net a labelled graph and an accompanying set of node matrices 
that satisfy Eqs.(2.1), (2.2) and (2.3), but don't necessarily satisfy the overall nor- 
malization condition Eq.(2.4). An acyclic CB pre-net always satisfies Eq.(2.4), and a 
paracyclic CB pre-net may not satisfy Eq.(2.4). This is why Following the procedure 
just discussed, the nodes of any acyclic pre-net can be relabelled chronologically as 
x 1; x 2 , • • • , xn_. The relabelled graph, even if it is not fully connected, corresponds 
to the right hand side of Eq.(2.13), except that some of the conditional probabili- 
ties on the right hand side of Eq.(2.13) might include redundant conditioning. (If 
a conditional probability P(x\y,R) is known to satisfy P(x\y,R) = P(x\y), then we 
would say that the expression P(x\y,R) shows redundant conditioning on R.) This 
correspondence between any acyclic pre-net and the right hand side of Eq.(2.13) guar- 
antees that acyclic pre-nets will satisfy the overall normalization condition Eq.(2.4). 
On the other hand, as we've shown with the 3-node pre-net of Fig.4, an iV-node 
paracyclic pre-net may not reduce to unity upon summing over its free indices. If one 
tries the process of "adding over all the current external arrows" on a pre-net with 
a cycle embedded in it, at some point the process comes to a stop and cannot be 
completed due to a lack of a current external arrow to sum over next. Note that, in 
some sense, Eq.(2.13) embodies the principle of causality. Thus, acyclic CB pre-nets 
preserve causality, and paracyclic CB pre-nets violate this principle. If one considers 
only acyclic graphs, as we shall do henceforth, then there is no difference between CB 
nets and CB pre-nets. 

Note that if one sums both sides of Eq.(2.11) over y and divides by P(x), one 
obtains 

P(z\x) = '£P(z\y)P{y\x). (2.14) 

y 

This last equation, vàlid for a Markov chain (x,y,z), is the so called Chapman- 
Kolgomorov equation. Equation (2.14) is represented diagramatically in Fig.7. If one 
sums both sides of Eq.(2.12) over y and divides by P(x) ) one obtains 
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P{z\x)= J £P(z\y,x)P(y\x). (2.15) 
y 

The last equation is a gener alization of the Chapman-Kolgomorov equation to arbi- 
trary random variables x,y,z triat don't necessarily form a Markov chain. Equation 
(2.15) is represented diagramatically in Fig.8. Note that in both Figs.7 and 8, we 
follow a process of adding over some of the arrows of a net to obtain a new net with 
fewer nodes. This process may be called coarsening or data compression, because it 
reduces the number of nodes and because the joint probability distribution of the new 
net carries less information than the joint probability distribution of the old one. In 
general, if we start with iV nodes x x , x 2 , ■ ■ ■ , x N , and we sum over xi, x 2 , ■ ■ ■ , x k , then 
the resulting probability distribution of the variables Xk+i,Xk+2, • • • ,%n can always 
be represented by a fully connected net with nodes x k+ i,x k+2 , ■ ■ ■ ,x N . 

In modelling a classical physical situation by a CB net, if one knows very 
little about the nature of the probability distributions involved, one can always use 
a fully connected net. Later on, if one learns that certain pairs of random variables 
are conditionally independent, one may remove the arrows connecting these pairs of 
variables without changing the value of the full net. 

Once we have designed a net architecture and the net's node matrices have 
been calculated, what next? How to use this information? The information is suffi- 
cient for calculating the joint probability P(xi, x 2 , • • ■ , xn), where 
the nodes of the net. From this joint probability, one can calculate P(x a \xb, x c , ■ ■ •), 
the probability that any one node x a of the net assumes one of its states x a , condi- 
tioned on several other nodes Xj,, x c , • • ■ of the net assuming respective states Xb, x c ,---. 
To go from P(x\, x 2 , ■ ■ ■ , xjy) to P(x a \xb, x c , • • ■), one can sum P(xi, x 2 , ■ • • , x^) over 
those variables which are not in the list x a , Xb, x c , ■ ■ ■ to get P(x a , Xf,,x c ,---). One 
can obtain P(xb, x c , ■ ■ •) similarly, and then divide the two. However, this brute force 
method takes no advantage of the particular topology of the net, and so it is very 
labor intensive. Artificial intelligence and expert systems workers have found an al- 
gorithm for calculating P(x a \xb, x c , ■ ■ •) from P(x\, x 2 , ■ ■ • , xn) that takes advantage 
of the net topology to reduce dramatically the number of arithmetical operations re- 
quired. For a discussion of this great achievement, the reader can consult Refs.|ÏJ-[||. 
The software package of Ref . implements this fast algorithm. 

3. EXAMPLES OF CB NETS 

The references at the end of this paper present examples of CB nets used in medical 
diagnosis @],@, monitoring of processes]]]], genètics 0, etc. And this is just a small 
fraction of the possible applications of CB nets. Indeed, any probabilistic model can 
be discussed in terms of CB nets. In this section, we do not intend to present a com- 
prehensive collection of CB net examples, but only to give a few simple introductory 
ones. 
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(a) DIGITAL CIRCUITS 

Figure 9 shows a Bayesian net version of an AND gate. The random variables 
x, y and z are binary (i.e. they assume vàlues in {0, 1}.) The node matrices for this 
net consist of the prior probabilities P(x), P(y) and of the conditional probability 
P(z\x,y). P(z\x,y) is given by the table in Fig.9. With states and 1 standing for 
false and true, respectively, this table is what one would expect if z — ( x and y ). 
Note that the columns of the node matrix P(z\x,y) add up to one, as required by 
Eq.(2.2). A node matrix whose entries are all either or 1 (like P(z\x,y) in this 
example) will be called a deterministic node matrix. 

CB nets representing OR and NOT gates can be defined analogously to the 
AND net. Then one can construct CB nets for any combination of AND, OR and 
NOT gates. 

(b) CONSTRAINT NODES 

Figure 10 shows a net for the binary random variables x, y and the random 
variable z G {0, 1, 2}. To fully specify this net, one must give the prior probabilities 
P(x), P(y) and the conditional probability P(z\x,y). The node matrix P(z\x,y) is 
given by the table of Fig.10. This node matrix is deterministic and agrees with 
z — x + y. Note that if we know for sure that x + y = 1, then we can "fix" the 
sum node z to one, and calculate probabilities, like P{x\z_ = 1), which have z = 1 as 
evidence. In such a case, we call z a constraint node, because it is used to enforce the 
constraint x + y = 1. 

As another example of a net that possesses a node which is useful to constrain, 
consider Fig.ll. In this net the random variables x, y and z are binary. To fully specify 
the net, one must give prior probabilities P{x), P(y) and the conditional probability 
P(z\x, y). P(z\x,y) is given by the table in Fig.ll. This table is what one would 
expect if z = (if x then y) . Since an "if x then y" statement does not say anything 
about what to do when x = 0, the probabilities P(z\x,y) are arbitrary and not 
necessarily deterministic when x = 0. The net of Fig.ll could be used by fixing the 
if-then node z to one and considering only probabilities with z = 1 as evidence. In 
such a case, we would say that the z node was a constraint node. 

(c) CL AUSER-HORNE EXPERIMENT 

Consider the Clauser-Horne experiment, which is used to observe violations 



of Bell-type inequalities|ï2|. Two partides, call them 1 and 2, are created at a 
common vèrtex and fly apart. Let À represent the "hidden variables". Particle 1 is 
subjected to a spin measurement along the direction A with outcome xf G {+, — } 
or along the direction A' with outcome xf G {+, — }. Particle 2 is subjected to a 
spin measurement along the direction B with outcome G {+, — } or along the 
direction B' with outcome G {+, — }. One can draw the net of Fig.12, which has 
nodes À, x_i and x^ . (Figure 12 really represents 4 nets, one for each of the following 
possibilities: (0 1} 2 ) = (A, B), (A, B'), (A', B), (A', B')). To specify the net, one must 
give probabilities -P(A), P(x^ 1 |À) and P(x2 2 |A). Note that if one sums the net of 
Fig.12 over A, one gets 
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P{x{\4>) = £P(^|A)P(^|A)P(A) , (3.1) 

A 

which is the starting point in the derivation of the Bell inequalities for the Clauser- 
Horne experiment. 

As a slightly more complicated experiment, one might choose at random 
whether 9\ = A or A' and whether 9 2 = B or B' . Such an experiment can be 
represented by a CB net with nodes 9_ l) x_±, 6_ 2 , x 2 and A, where Xj for j G {1,2} rep- 
resents the outcome of a measurement on particle j. See Fig.13. This net is specified 
if we give probabilities P(xi\9i, A), -P(#i), P(x 2 \9 2: A), P{9 2 ) and -P(A). Summing the 
net over A and dividing by P(9i)P{9 2 ) yields an equation analogous to Eq.(3.1): 

p(x 1 ,x 2 \e 1 ,e 2 ) = Y / P{x 1 \e 1 ,x)P(x 2 \e 2 ,x)p(\) . (3.2) 

A 

(d)RANDOM WALK 

Suppose that a particle moves in a straight line, taking unit length steps either 
forwards or backwards, with probabilities p + and respectively, where p + +p_ = 1. 
Let Xj be the position of the particle at time j, with j G {0, 1, 2, • • • , }. Assume Xo, 
the starting position, is zero. Define Axj by Axj = Xj — Xj-i, for j e {1, 2, • • • , }. 
Figure 14a shows a CB net that represents the probability distribution of Xq,^, ■ ■ ■ 
and Ax^Ax-,. ■ ■ •■ Clearly, ^ G {0, ±1, ±2, • • • , ±j} and Axj G {±1}. To fully 
specify the net of Fig.l4a, we must give the probability matrices associated with each 
node. These matrices are 



P(x = 0) = 1 , 

P{Axj±l)=P± forj = l,2,... , (3.3) 
P(x.j = ylxj-! = x, Axj ±l) = 6(y,x±l) for j = 1, 2, • • • , 

where 5 is the Kronecker delta function. By summing the net of Fig.l4a over all 
possible vàlues of the indices Axj that label the arrows from the Axj to the x_j nodes, 
one obtains Fig.l4b. For this coarser net, one obtains 

= 0) = 1 , (3 4) 

P{x.j = = x) = p + 5{y, x + 1) + pS(y, x - 1) for j = 1, 2, • • • . 

One may go even further: For any k G {0, 1, 2 • • •}, one may sum Fig.l4b over all Xj 
such that j £ {0, k}. One obtains the net of Fig.l4c, with 

P(x = 0) = 1 , (3.5a) 
P(x k = Xk \x = 0)=( r + S ) f + p s _ , (3.5b) 
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where r — s = Xk and r + s = k, and the first factor on the right hand side of Eq. (3.5b) 
is a combinatorial factor. 



4. THEORY OF QB NETS 

In this section, we will define QB nets and explain how to use them to calculate 
quantum mechanical conditional probabilities. The next section will give examples 
of QB nets. 

Like a CB net, a QB net consists of two parts: a labelled graph and a collec- 
tion of node matrices. These two parts must satisfy certain constraints that will be 
specified below. 

External and infernal arrows, external and infernal nodes, and root nodes are 
all defined in the same way for QB nets as for CB nets. All nodes of a QB net must 
be either infernal or external. 

We define each node of a QB net to represent a numerical value, and the 
whole net to represent the product of all these node vàlues. We assign a numer- 
ical value to each node as follows. First, we assign a random variable to each 
node. Suppose the random variables assigned to the N nodes are x_ ^ , X o , • • • , x m . 
Define Zn = {1, 2, • • • , N}. For any finite set S, let \S\ be the number of ele- 
ments in S. If S — {ki, Jï2, ■ • • , k\s\} C Zn, define (x.)s = (xku %k 2 , • • • , x k lsl ) and 
{x-)s — (xj ei ,Xj e2 ,···,x_i e . s .). Sometimes, we also abbreviate (x.)z N (i-e., the vector 
that includes all the possible Xj components) by just x., and (x.)z N by just x. . For 
j G Zn, we imagine node x_j to lie in state Xj. We also imagine all arrows leaving the 
node Xj to lie in state Xj, and thus we label all of them Xj. At this point we've shown 
how to label each arrow in the graph by x^ for some k G Z^. Define Sj to be the 
set of all k such that an arrow labelled x^ enters node x_j. Now we assign a complex 
number .<4[xj|(2;.)s.] to node Xj. i4[xj|(x.)s.] is the node matrix for node xf Xj is 
the matrix's row index and (x^Sj is its column index. We require that the quantities 
A[x_j|(x.)s.] be probability amplitudes that satisfy 

Y^\A[x 3 \{x.) Sj ]\ 2 = l : (4.1) 

Xj 

where the sum in Eq.(4.1) is over all states Xj that the random variable x_j can assume, 
and where Eq.(4.1) must be satisfied for all j G Z N and for all possible vàlues of the 
vector (x.)s. of random variables. The QB net is taken to represent the product of all 
the probability amplitudes i4[xj|(x.)s.] for j G Z N . This product is a function A(x.) 
of the current states xi, x 2 , • ■ • , x^ of the nodes. Thus, 

A{x.)= I] A[ Xj \(x.) 8j ] . (4.2) 
jez N 

Let Z e ^ 1 be the set of all j G Z N such that Xj is an external node, and let Z™ 1 be 



11 



the set of those j G Z N such that Xj is an internal node. Clearly, Z™ 1 and Z 1 ^ 1 are 
disjoint and their union is Z N . We require A(x.) to satisfy 



E 

(x.) z ea 



E 

(a\ ) z int 



(4.3) 



and 

El^-)| 2 = 1- (4-4) 

X. 

Note that as a consequence of Eqs.(4.1) and (4.4), given any QB net, one can 
construct a special CB net by replacing the value -AfxjKx.),^] of each node by its 
magnitude squared. We call this special CB net the parent CB net of the QB net 
from which it was constructed. We call it so because, given a parent CB net, one can 
replace the value of each node by its square root times a phase factor. For a different 
choice of phase factors, one generates a different QB net. Thus, a parent CB net may 
be used to generate a whole family of QB nets. 

A QB pre-net is a labelled graph and an accompanying set of node matrices 
that satisfy Eqs.(4.1), (4.2) and (4.3), but don't necessarily satisfy Eq.(4.4). A QB 
pre-net that is acyclic satisfies Eq.(4.4), because its parent CB pre-net is acyclic and 
this implies that Eq.(4.4) is satisfied. If one considers only acyclic graphs, as we shall 
do henceforth, then there is no difference between QB nets and QB pre-nets. 

In the second quantized formulation of quantum mechanics, one speaks of M 
modes represented by M annihilation operators ai, a 2 , • • • , om, which obey certain 
commutation relations amongst themselves. Define Zq + = {0, 1,2, • • •} and let rii G 
Z 0+ for i — 1 to M. For M modes, one uses quantum states like 

_tm „t«2 n \n M 

«i «2 : | ). (4.5) 



\Jn M \ 

The state given by Eq.(4.5) is specified by a vector (ni, n 2 , • • • , n M ) of occupation 
numbers. For the rest of this paper, we will consider only QB nets whose states Xj 
are vectors (7ij,i, 7ij,2, • • • , n j,K ) of occupation numbers. Define T to be the set of all 
a such that n a is an occupation number of the QB net under consideration. For 
example, a = (j, 1) in n^i, where is a component of Xj = (rij^, n Ji2 , • • • , nj,Kj)- 
For any set T' = {«i,a 2 , • • • ,«|r'|} C T, let (n.) r = (n ai ,n a2 , ■ ■ ■ ,n Q|r/| ) and (n.) r , = 
(z^qj, 21a 2 , 21a, r /|)- We will sometimes abbreviate (n.)r by n. and (n.)r by n. . 
For each j, define Tj to be the set such that (^Or^. is the occupation number vector 
specifying the state of node Xj. Thus, Xj = {n.) rj and Xj = {n.)r r If 

Xj IS cl II external 

node, call all the components of (n.)r 3 external occupation numbers, and if Xj is an 
internal node, call all the components of (n.)r\, internal occupation numbers. Let Y ext 
be the set of all a such that n a is an external occupation number, and let T mt be the 
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set of all a such that n a is an internal occupation number. Clearly, T mt and T ext are 
disjoint and their union is V. Note that (x.) Z ext = (n.)r^*, (x.) Z int = (n.) r mt, and 
( x -)z N = ( n -)r (Analogous statements with x and n underlined to indicate random 
variables also hold). 

Consider a classical probability distribution P{n.). For any r C T, we may 
define the characteristic function Xc by 

Xc[(n.)r ] = E P K) II ó ( m «> "«) , (4.6) 
in. aer 

where 6 is the Kronecker delta function. (The "c" subscript in Xc stands for "clas- 
sical".) If T H and T E are disjoint subsets of T, one defines the classical conditional 
probability that (n.)r H = (n.)r H (hypothesis) given that (n.)r E = {n'.)r E (evidence), 
by 

P[(n.) rH = (n.) rH \(n.) TE = (n[) rE } = ■ (4.7) 

Xc[(n.)r E \ 

By setting T E = 4> in the last equation, we conclude that for any T C T, 

P[(n.)r = (n.) ro ] = Xc[(n.)r ] • (4-8) 
According to the definition Eq.(4.6), if T and T' are disjoint subsets of T, then 

E Xc[(n.)ro,(níV ]=Xc[(n.)r ]- (4-9) 
K.) r; , 

Equations (4.8) and (4.9) imply that 

E P[(n.)r ,(n'.h } = P[(n.)r ], (4.10) 

for T , Tg C T, T n T' = <fi. Equation (4.9), when applied to Eq.(4.7), yields 

E P[(nOrJ(ní)rJ = l. (4.11) 

(».)r H 

Furthermore, from Eq.(4.7) it is obvious that 

P[(n.) rH \(n'.) rE }>0. (4.12) 

Equations (4.11) and (4.12) imply that Eq.(4.7) is an adequate definition of condi- 
tional probability. 

Now consider a quantum mechanical probability amplitude A(n-). For any 
T C T, we define the characteristic function x by 
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Xl(n.)r }= E E A(m.) H 5(m a ,n a ) . (4.13) 

(m.) r exí (m-) r *™í <*er 

If Th and are disjoint subsets of T, one defines the quantum mechanical conditional 
probability that (r.)y h = (n.)r H (hypothesis) given that (n.)r E = ( n '-)r E (evidence), 
by 



P[(n.) rH = (n.) TH \(n.) rE = (n'.) rE ] = ^^f^L i - (4-14) 

Note that the denominator of the right hand side of Eq.(4.14) depends on both T E and 
IV, unlike the analogous denominator in the classical definition Eq.(4.7). In quantum 
mechanics, with \c replaced by Xi Eq.(4.8) is not true for all T C T (although it 
is true if T C Y ext ). In quantum mechanics, with \ c replaced by Xi Eqs.(4.9) and 
(4.10) are not necessarily true, whereas Eqs.(4.11) and (4.12) are obviously always 
true. Because Eqs.(4.11) and (4.12) are satisfied in quantum mechanics, Eq.(4.14) is 
an adequate definition of conditional probability. 

Note that both the classical and quantum mechanical definitions of conditional 
probability Eqs. (4.7) and (4.14) implicitly assume that we perform non-destructive 
measurements on any infernal nodes that might be measured. Hence, if a particle 
is detected at an infernal node, it is allowed to continue past that node so that it 
can reach other nodes further downstream. If one wishes to perform a destructive 
measurement on a particle when it passes through an infernal node x, then one is 
representing the physical situation by the wrong CB or QB net; what is required is a 
net that has the node x as an external node. 

In Appendix A, the classical and quantum mechanical definitions of conditional 
probability Eqs. (4. 7) and (4.14) are generalized so that they allow either sharp or 
fuzzy hypotheses and pieces of evidence. If the evidence narrows the set of possible 
vàlues for n x to a single number, say n\ — 1, then we say that the evidence for n x is 
sharp. If the evidence doesn't do this, then we say that the evidence for n x is fuzzy. 
Sharp and fuzzy hypotheses are defined analogously. 

Note that if we had used Eq.(4.7) (with Xc replaced by x) as the quantum 
mechanical definition of conditional probability, there would have been no guarantee 
that Eq.(4.11) would be satisfied. For given T H ,T E and (n'.)r E , define the quantum 
non-additivity factor f qna by 

t rr r < >\ 1 ^(™-)r H x[(™.)r H , (n[) rE ] (A , ^ 

f q na[T H , T E , (n.)r E \ = \\ (rt ,\ ï ■ ( 4 - 15 ) 

This quantity will be calculated for the examples in the next section. If f qna = 1, then 
Eq.(4.7) (with Xc replaced by x) an d Eq.(4.14) agree; if f qna ^ 1, then Eq.(4.7) does 
not give a well defined probability distribution for (ji.) Th whereas Eq.(4.14) does. 
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In quantum mechanics, x[( n -)r H , ( n '-)r E ] is proportional to the number of oc- 
currences of (n-)r H = { n -)r H and (n.)r E = {n'.)r E in an experiment that measures 
all the Th and Te nodes. x[(w'.)r B ] is proportional to the number of occurrences of 
(ïL-)r E = { n '-)v E in an experiment that measures only the Y E nodes, leaving the Y H 
nodes undisturbed. Thus, Eq.(4.7) (with Xc replaced by x) refers to the ratio of the 
number of occurrences in two different types of experiments (one type measuring the 
ThUTe nodes and the other only the Te nodes). On the other hand, Eq.(4.14) refers 
to the ratio of occurrences in a single type of experiment (measuring the U Te 
nodes). It seems to us that the latter ratio is the more useful of the two. 

By using CB and QB nets, one is led easily and naturally to calculate proba- 
bilities by considering sums over paths, the type of sums advocated by Feynman for 
quantum mechanics and by Kac for Brownian motion. Indeed, one can express the 
classical and quantum mechanical definitions of conditional probability Eqs.(4.7) and 
(4.14) in terms of sums over paths rather than sums over node states. We do so for 
arbitrary CB and QB nets in Appendix B. Specific examples illustrating the afhnity 
of QB nets with sums over paths can be found in the next section and in Appendix 
C. Appendix C presents a QB net which yields the Feynman path integral for a single 
mass particle under the influence of an arbitrary potential. Using an approach similar 
to that of Appendix C, it should be possible to define QB nets that yield the Feynman 
integrals employed in non-relativistic and relativistic quantum field theories. 

5. EXAMPLES OF QB NETS 

In this section, we will present the results of a computer program that calculates 
conditional probabilities for QB nets. In particular, we shall consider QB nets for 
experiments containing either two or three Stern-Gerlach magnets. We will restrict 
our attention to experiments involving a single particle. For simplicity, for each 
experimental configuration, we will assume that the magnètic field vectors of all the 
Stern-Gerlach magnets are coplanar. 

Figures 15, 16 and 17 show the three kinds of nodes that will be used in this 
section. 

In Fig.15, the triangle represents a root node. This node will stand for ip nz _ nz+ , 
with (n z _,n z+ ) G {(0, 1), (1, 0)} and |^ i| 2 + |^io| 2 = 1- ^n z _n z+ is just the initial 
wavefunction for the single particle under consideration. 

In Fig.l6a, the black-filled circle represents a marginalizer node. The single 
incoming arrow is in a state characterized by a vector (jii, ri2, • • • , nx) of occupation 
numbers. The outgoing arrow is in a state characterized by a single occupation 
number n\. The amplitude associated with the node is 

A(n' 1 \n 1 ,n 2 , ■ ■ ■ ,%) = S(n' 1: ni) . (5.1a) 
Thus, a marginalizer node takes a vector of occupation numbers and projects out one 
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of its components. Note that Eq.(5.1a) satisfies Eq.(4.1). 

In Fig.l6b, the black-filled circle with a phase factor next to it represents a 
phase shifter node. The single incoming (outgoing) arrow is in a state characterized 
by a single occupation number n\ (n[). The amplitude associated with the node is 

Airí^ri!) = e*5(ní,ni) , (5.1b) 

where £ is a real constant. Note that Eq.(5.1b) satisfies Eq.(4.1). 

In Figs.17, the white-filled circles represent Stern-Gerlach magnets. The out- 
going arrows are labelled by a vector n™ = (n"_,n" + ) of occupation numbers. Here 
the unit vector ü represents the direction of the magnet's magnètic field, and a la- 
bels the magnet ( ü is not enough to label the magnet if the experiment contains 
more that one magnet whose magnètic field points along the ü direction.) The vector 
n" G {(0, 0), (0, 1), (1, 0)} specifies a state (see Appendix D) 

«!)<- (a&rl·lO) • (5-2) 
The creation operators a"í_ and a„+ create partides in the |— u ) and \+ u ) states 
respectively. In Fig.l7a there is one arrow entering the node, whereas in the Fig.l7b 
there are two. In general, since we are considering a single particle experiment, there 
may be any number of arrows entering the node, but only one may be in state 1, 
all others must be in state 0. In Figs.l7a and 17b, the amplitude assigned to each 
node is given by a table next to the graph. Clearly, the tables in Figs.l7a and 17b 
both satisfy Eq.(4.1). For each experimental configurat ion, we will assume that the 
magnètic fields of all magnets are coplanar. The plane containing these magnètic field 
vectors may be chosen to be the X-Z plane with = 0. Thus, the matrix elements in 
the tables of Figs.l7a and 17b are given in terms of angular parameters by Eqs.(D.7) 
and (D.8). 

(a) EXPERIMENTS WITH 2 STERN-GERLACH MAGNETS 

We will consider 2 configurations with 2 magnets: Fig.18 (the tree diagram) 
and Fig.19 (the loop diagram). The diagrams are like road maps, with the arrows 
representing the various roads along which the particle may travel. 

A single particle travelling through the experimental configuration of Fig.18 
could exit through either the n z _, n u _ or n u+ nodes. There is only one possible path 
leading to each of these outcomes. Thus, one has 

Fh = A(tti) , 

A(m) = V>io , 

Fh = A(n 2 ) , 

A(n 2 ) = (- u \+ g )ip 01 , 
FI 3 = A(tt 3 ) , 

A(n 3 ) = {+ u \+ z )^oi ■ 

For i G {1,2,3}, A(iTi) is the amplitude for path 7Tj. For j G {1,2,3}, FIj is a 
Feynman integral; that is, the sum of the amplitudes for all paths with a given final 
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(5.3) 



state. (See Appendix B). We have already checked that Eq.(4.1) is satisfied. We did 
so when we defined the wavefunction, marginalizer and Stern-Gerlach nodes. It is 
easy to show that Y^=i \FIj\ 2 = 1, so Eq.(4.3) is also satisfied for this net. 

A single particle travelling through the experimental configuration of Fig.19 
could exit through either the n u _ or the n u+ nodes. There are two possible paths 
leading to each of these two outcomes. Thus, one has 



Fh = A(tti) 


+ A(n 2 ) , 


AM = <- 


u 


-z)lj>10 


A(n 2 ) = (■ 


~u 




FI 2 = A(n 3 ) 


+ A{tc 4 ) , 


M*s) = (- 


fui 




A(n 4 ) = (- 


fui 


+z)Í>0l 



(5.4) 



As before, for i e {1, 2, 3}, -A(7Tj) is the amplitude for path 7Tj, and for j 6 {1, 2}, F/j 
is a Feynman integral. It is easy to show that Y^j=i \FIj\ 2 = 1, so Eq.(4.3) is satisfied 
by this net. 

We will call an arrow or a node simple if its state is characterized by a single 
occupation number (for example, a marginalizer node and its outgoing arrow are both 
simple). We will say that a net is fully marginalized if the incoming arrows of any 
node that is not a marginalizer are all simple, and all external arrows are simple. The 
nets of Figs.18 and 19, and, in fact, all the nets considered in this section, are fully 
marginalized. 

For each of the nets of Figs.18 and 19, we used a computer program to calculate 
probabilities with one or two hypotheses and with zero, one or two pieces of evidence. 
More precisely, for the 2 QB nets in Figs. 18 and 19, and for their parent CB nets, for 
n, n', m and m' G {n u+ ,n u _,n z+ ,n z _}, for all n, n', m and w! G {0, 1}, we calculated 
(using Eqs.(4.7) or (4.14)) the unconditional probabilities P(n) and P(n,n'), and the 
following conditional probabilities: 

P(n = n\m = m) , (5.5) 
P(n = n\m = m,rrí = m!) , (5.6) 
P(n = n,rí = n'\m = m), (5.7) 

P(n = n,rh=n'\m = m,rn/ = m'). (5.8) 

The table of Fig.20, call it the evidence- case file, gives the various sets of 
evidence that were considered. For example, in case 2 (the row that starts with a 
2), we assumed n z+ = 0, whereas the vàlues of the remaining occupation numbers ( 
) were assumed to be unknown (i.e., we assumed there was no evidence 
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as to whether they were or 1 and this is indicated in Fig.20 with a blank space ). 
In case 11, we assumed n z _ = 1 and n z+ = 0, whereas the vàlues of the remaining 
occupation numbers ( n u+ ,n u _) were assumed to be unknown. Note that for the 
evidence-cases 2 to 9 we considered 1 piece of evidence (as in Eqs.(5.5) and (5.7)), 
and for the evidence-cases 10 to 33 we considered 2 pieces of evidence (as in Eqs.(5.6) 
and (5.8)). 

To get numerical vàlues for the probabilities associated with the nets Figs.18 
and 19, particular vàlues had to be assumed for the initial wavefunction and for the 
magnètic field direction of each magnet. We used 

1 + Í 1 7T 

^oi = ~y~ , ifoo = -j= , Q U -0 Z = - . (5.9) 

For the tree graph of Fig.18, all the probability distributions for the QB net 
and for its parent CB net were identical. This is not surprising since tree graphs, i.e., 
graphs without any loops, have no interfering paths. That is, for a given final state, 
there is only one possible path that produces that final state. More interesting results 
were obtained for the loop graph of Fig.19. 

Figures 21 and 22 show our computer program's output for the loop net, for 
cases 1, 2, 4, 10, 12 (See evidence-case file Fig.20). Figure 21 gives one-hypothesis 
probabilities of the type Eq.(5.5) and (5.6), whereas Fig.22 gives two-hypotheses 
probabilities of the type Eq.(5.7) and (5.8). 

First consider Fig.21. Columns A to D refer to the parent CB net and columns 
F to I to the QB net. Column A for the CB net (F for the QB net) gives the identity 
of the occupation number n in Eqs.(5.5) and (5.6). Columns B and C for the CB net 
(G and H for the QB net) give the probabilities that n = and n = 1, respectively, 
in light of the evidence. Column D for the CB net (I for the QB net) gives the 
quantity f qna defined by Eq.(4.15). Note that for evidence-case 10, there was no 
output, because the computer program detected a contradiction. (In case 10, we 
were assuming that n z+ = n z _ = 0, which is impossible since n z+ + n z _ = 1.). More 
interesting is the evidence-case 4. In this case, for the QB net, column I indicates 
that f qna ^ 1 for the distributions P{n z+ \n u+ = 0) and P{n z _\n u+ = 0).[|I3| 

Next consider Fig.22. Columns A to F refer to the parent CB net and columns 
H to M to the QB net. Column A for the CB net (H for the QB net) gives the 
identity of the occupation numbers n and n' in Eqs.(5.7) and (5.8). Columns B 
to E for the CB net (I to L for the QB net) give the probabilities that (n,n') = 
(0, 0), (0, 1), (1, 0), (1, 1), respectively. Column F for the CB net (M for the QB net) 
gives f qna . For example, in Fig.22 we see that for evidence-case 4, for the QB net, all 
the probability distributions except P(n u+ , n u _\n u+ = 0) have f qna ^ 1. JHf 

(b) EXPERIMENTS WITH THREE STERN-GERLACH MAGNETS 

We will consider 7 configurat ions with 3 Stern-Gerlach magnets: Figs.23 to 

29. 

For Fig.23, define 
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Fh = A(m) , 
A(m) = íp 10 , 

FI 2 = A(7T 2 ) , 

A(n 2 ) = (- v \+ z )ip 01 , 
Fh = A(ti 3 ) , 

A M = (-u\+v) ( + v\+z)^01 , 

FI 4 = A{n A ) , 

A(lT 4 ) = (+ u \+v)(+v\+z)^01 ■ 

It is easy to show that X^=i \FIj\ 2 = 1, so Eq.(4.3) is satisfied by this net. 
For Fig.24, define 

Fh = A(m) + A(rr 2 ) , 

A M = (-v\-z)ipio , 

A(n 2 ) = (- v \+z)^oi , 
FI 2 = A{n 3 ) + A{tt 4 ) , 

A M = (-u\+v)(+v\-z)^W , 

A(tt 4 ) = (- u \+ v ){+ v \+ z )ipoi , 
FI 3 = A(n 5 ) + A(tt 6 ) , 

A M = (+v.\+v)(+v\-z)l/>10 , 
A M = (+v.\+v)(+v\+z)^01 ■ 

It is easy to show that Y%=i \FIj\ 2 = 1, so Eq.(4.3) is satisfied by this net. 
For Fig.25, define 

Fh = A(tti) , 
A(m) = íp w , 

FI 2 = A{n 2 ) + A{n 3 ) , 

A(iï 2 ) = (- u |-„)(-„|+ z )V'oi , 

M^) = (-u\+v) (+v\ +z)ïp01 , 

FI 3 = A(ir A ) + A(ir 5 ) , 

A(7T 4 ) = (+ u |-t;)(-t;|+z)^Òl , 
A(7T 5 ) = (+ u \+v)(+v\+z)thl ■ 

It is easy to show that Y%=i \FIj\ 2 = 1) so Eq.(4.3) is satisfied by this net. 
For Fig.26, define 
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(5.14) 



Fh = Aim) + A{n 2 ) + A{ir 3 ) + A{n 4 ) , 

^("l) = (-u\+v) (+v\-z)lplO , 
A{lï 2 ) = (-u\-v)(-v\-z)Ao , 
A M = (-u\+v)(+v\+z)ïp01 , 

A(n 4 ) = (- u \-v){-v\+z)lp01 , ( r -,oN 

FI 2 = A(n 5 ) + A(rv 6 ) + A(tt 7 ) + A(tt s ) , {ÚAà) 

Afa) = (+u\+v)(+v\-z)ïpW , 
A M = (+u|-t;)<-t;|-«)^10 , 

A(n 7 ) = (+ u \+v)(+v\+z)ïpoi , 
A(ir 8 ) = (+ u \-v)(-v\+z)4>oi ■ 

It is easy to show that J2j=i \FIj\ 2 = 1, so Eq.(4.3) is satisfied by this net. 
For Fig.27, define 

Fh = Afa) , 

A M = (-v\-z)lplO , 

FI 2 = A(n 2 ) , 

A(ir 2 ) = (+ v \- z )ip 10 , 
FI 3 = A(tt 3 ) , 

A M = (-u\+z)^01 , 

Fh = A{it 4 ) , 

A(n 4 ) = (+ u |+ z )^ i • 

It is easy to show that Yfj=\ Wh\ 2 = 1, so Eq.(4.3) is satisfied by this net. 
For Fig.28, define 

Fh = A(tti) , 

A ( n l) = (-v\-z)lpl0 , 

FI 2 = A{n 2 ) + A{n 3 ) , 

A(n 2 ) = e^(+ u \+ v )(+ v \- z )ÍJ 10 , 

^(^s) = (+«1+2)^01 , 

FI 3 = A(ir 4 ) +A(ir 5 ) , 

A(tt 4 ) = e^(- u \+ v )(+ v \- z )i) W , 

= (-u|+«)^01 • 

Note that we have included a phase factor which can be thought of as arising from 
the Stern-Gerlach node u when its input n v+ equals 1 (or from node v when its output 
n v+ equals 1, or from node z when its output n 2 _ equals 1). Alternatively, the phase 
factor may be though of as arising from a phase shifter node, not pictured in Fig.28, 
located in the middle of the arrow labelled n v+ . If £ = 0, then Y%=i \FIj\ 2 1 an d 
therefore Eq.(4.3) is violated. On the other hand, if 

e* = , (5.16) 



(5.15) 
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then £j=i \FI Ó \ 2 = l.[[T§ 
For Fig.29, define 



FJx = A(7n) + A(tt 2 ) + A(tt 3 ) , 

^(^l) = (+u|+*)^01 , 

A(n 2 ) =e i t(+ u \+ v )(+ v \- z )ÍJ W , 

A(vr 3 ) = e*(+ u \- v )(- v \- z )i/) W , . , 

F/ 2 = A(tt 4 ) + A(vr 5 ) + A(tt 6 ) , 

A(7T 4 ) = (- U |+ 2 )^ 1 , 

A(tt 5 ) = e^(- u \ +«)(+„!-*) V»io , 



A(7r 6 ) = e ïÇ (- u |- í ,)(- í ,|- 2 )^ 
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It is easy to show that Z)j=i |F/j| = 1 for any value of £. However, we believe the 
correct choice of £ to be the one given by Eq.(5.16). For this choice, P{n u _ = l\n v _ = 
0) is the same for Figs.28 and 29. And this is what we expect. Otherwise, empty 
De Broglie waves could influence the outcome of an experiment (and thus could be 
detected), which does not appear to be the case experiment ally. 

For the QB nets of Figs.23 to 29, and for their parent CB nets, for n, n', m and 
ïïl' G {n u± , n v ±, n z± } , for all n, n', m and m! e {0, 1}, we calculated the conditional 
probabilities of Eqs.(5.5) to (5.8). 

The table of Fig.30 is an evidence-case file that gives the sets of evidence that 
were considered for the nets of Figs.23 to 29. For example, in case 2 ( the row that 
starts with a 2), we assumed n z+ = 0, whereas the vàlues of the remaining occupation 
numbers (n u ±, n v ±, n z J) were assumed to be unknown. 

To get numerical vàlues for the probabilities associated with the nets of Figs.23 
to 29, particular vàlues for ^ 01 , ip 10 , 9 U , 9 V , 6 Z were assumed. 

For the tree graphs Figs.23 and 27, and also for the non-tree graph Fig.28, we 
found that the QB net and its parent CB net yielded identical probability distribu- 
tions. 



APPENDIX A. CONDITIONAL PROBABILITIES 
FOR FUZZY MEASUREMENTS 

In this appendix, we will generalize the classical and quantum mechanical definitions 
of conditional probability Eqs.(4.7) and (4.14) so as to include either sharp or fuzzy 
hypotheses and pieces of evidence. 

The following simple results from set theory are relevant. Given any finite set 
S, we will denote the number of elements in S by |ií>|. Given two sets R and S, we 
define the direct product set R x S by 

R x S = {(x,y)\x G R,y G S} . (A.l) 
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R x S is also denoted by the vector of sets (R, S). Given sets R±, R 2 , Si, S 2 , it is easy 
to show that (see Fig.A.l) 



(R 1 x R 2 ) n x S 2 ) = n Si) x (R 2 n S 2 ) . (A.2) 

Therefore, (Ri x i? 2 ) H (Si x S 2 ) = if and only if R± fi Si = or i? 2 H S 2 = 0. Given 
any set S and any integer n, we will denote by S n the product 5xSx--·xSofn 
copies of S. 

Suppose Q a C Z 0+ for each a. If T' = {ai,a 2 , • • • ,a\r>\} C T, we define 
the vector of sets or direct product set (Q-)r' by (Q-)r' = (Q ai ,Qa 2 r " ' ,Qa, r ,X 
Equivalently, one may write (Q-)r' = Q ai x Qa 2 x ' ' ' x Qa, r ,,- Note that (Q-)r' C Zq + '. 
Sometimes, we will abbreviate (<5-)r by just Q.. Two direct product sets (-R.W'i and 
(S.)|r'| are disjoint if and only if there exists an a G V such that i? a fi S a = <fi. 

We will need to consider collections H = {H!\i = 1,2 - ■ ■ ,\H\} of direct product sets 

Iri in 
H! C Zq+. Such a collection will be said to be a partition of Zo+ if any pair of distinct 

sets HÍ and H j . is disjoint and u|!Jifí = 4+- 

We are interested in P(n. e i/. |n. E E.), i.e., the probability that r?^ e ií a for 
all a E r given that np E -E/3 for all f3 E T. iJ Q C Zo+ is the hypothesis for n^, and 
-E/3 C Zo+ is the evidence fornp. If -£f Q contains only one element of Z 0+ , then we say 
that the hypothesis for n a is sharp, whereas if H a contains more that one element of 
Z 0+ , we say that the hypothesis for n a is fuzzy. Analogously, the number of elements 
in Ep determines whether the evidence for na is sharp or fuzzy. 

For any set S, define the indicator function ls(%) by 

1 , v f if x g S , 
!*(*) = | 1 [íxeS ■ (A.3) 

Consider first a classical probability distribution P(n.). One defines 



To write the last equation more succinctly, it is convenient to define the filter function 

iri 

/Q.(n-), for any direct product set Q- C Zq+, by 

/«>•) = n^w- ( A - 5 ) 

It is also convenient to define the characteristic functional Xc[K] of any function K{n.) 
of n- by 

X c [ff] = £P(n.)tf(n.). ( A -6) 
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(The c subscript in Xc stands for "classical"). Now Eq.(A.4) can be written succinctly 

as 



P(n. G H.\n. G E.) 



Xclfl·l. 



nE. 



Xc[f. 



E. 



(A.7) 



Note that fn.nE. — Íh.Íe.- Note that if H = {H!\i = 1, 2, • • • , \H\} is a partition of 
4+, then Eq.(A.7) satisfies 



\H\ 



J2 p (n- e Hí\n. G E.) = 1 . 



(A.8) 



Now consider a quantum mechanical probability amplitude A(n.). We define 
the characteristic functional x[K] by 



X[K}= £ 



E A(n.)*(n.) 



(A.9) 



Suppose 7í = {Hl\i = 1, 2, • • • , |7i|} is a partition of Z +. In quantum mechanics, the 
conditional probability that n. G H] given that n. & E., depends on the choice of 
partition 7í. We define this probability by 



P H (n. G Hí\n. G E.) 



nE. 



EjJi x[fH j r\E.1 



(A.10) 



Equation (A.8) is trivially satisfied by Eq.(A.lO). 



APPENDIX B. CONDITIONAL PROBABILITIES 
EXPRESSED IN TERMS OF PATH SUMS 

In this appendix, we will express the classical and quantum mechanical definitions 
Eqs.(A.7) and (A.10) of conditional probability in terms of sums over paths. Simple 
examples of the results of this appendix may be found in Section 5 of this paper. 

Let P(n.) and A(n.) be the vàlues of a CB net and a QB net, respectively. 
Suppose that V(n.) = P{n.) if a CB net is being considered and V(n.) = A(n.) if a 
QB net is. 

Let II be the set of all possible paths. II is defined so that there is exactly one 
7i G II for each n. that satisfies V(n.) ^ 0. Thus, there is a one to one onto map 
n.(7r) from II into {n. \V(n.) ^ 0}. 

For any V C T, if m. is such that (m.)r' = (n.)r>, call m. an extension of 
(n.)r'. Let £ be the set of all possible final states. S is defined so that there is 
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exactly one o G E for each (n.)p^t for which there exists an extension n. such that 
V(n.) ^ 0. Thus, there is a one to one onto map from E into {(n.)r^t\V(n.) ^ 0}. 

Define the partitíon function Z to be a function from II to E such that Z(ir) = a 
if the path tc has o as final state. For each o G E, let C(cx) be the set of all paths 
7T that have o as final state; i.e., C(a) = {n G IT|Z(7r) = a}. Clearly, the sets 
(equivalence classes) C(a) and C(cr') are disjoint when <r ^ a', and U CTe £C(cr) = II. 

For each n G II, define V(n) by 'P(tt) = V(n.(n)). 

For any direct product set Q- C Z +, define the filter function /q.(vt) by 

7 Q » = nw^wi- ( Ri ) 

«er 

Classically, one defines the characteristic functional xj-^l °f an y function K(tt) 

of 7T by 

Xct^l = E E ^(^(tt) . (B.2) 

Then 

P(n. G F. n. G £■) = _ . . (B.3) 

XcUe.) 

Quantum mechanically, one defines the characteristic functional ~x\K\ by 



X[K] = E 



o-es 



2 



(B.4) 



in 

If 7í = {H!\i = 1, 2, ■ • • , |7í|} is a partition of Z ^, then one defines 

xifl·lïnE. ) 



P n {n. G #?|n. G E.) = _ m{ _ : • (B.5) 



APPENDIX C. QB NET FOR 
SINGLE PARTICLE 
FEYNMAN PATH INTEGRAL 

In this appendix, we will present a QB net which yields the Feynman path integral |To| 
that in turn yields the Schroedinger equation for a single mass particle in an arbitrary 
potential. 

We begin by restricting positions x to lie in the interval ("box") [-ySf], 
where L is much larger that any other length (including the particle's position co- 
ordinate) occurring in the problem. Then we divide the box into subintervals of 
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infinitesimal length Ax, and call the midpoints of these subintervals x s with s G 
{0, ±1, ±2, • • • , ±N x /2} = Z x . Similarly, we restrict times t to lie in the interval 
[0, T\. Then we divide the interval [0, T] into subintervals of infinitesimal length At, 
and call the midpoints of these subintervals U with i G {0, 1, 2, • • • , N t } = Z t . 

Henceforth, for any function f(x s ), we will use /(•) to represent the vector 
whose sth component is f(x s ). For example, S Xr (-) will represent the vector whose 
sth component 5 Xr (x s ) is 1 if r = s and zero otherwise. 

As in Fig.C.l, we assign a node and a random variable n(x s ,ti) to all space- 
time lattice points (x s , U) with s G Z x and i G Z t — {0} and to the lattice point (x , t ). 
(Thus, at time í — 0, only position x gets a node). We draw arrows pointing from 
each node to all nodes occurring a time At later. We also draw external arrows 
pointing out of each node at time íjv t . The net starts at time to with a single root 
node at position xq, and it ends at time íjy 4 with external nodes at each position x s 
for all s G Z x . The random variables n(x s ,ti) are occupation numbers that assume 
vàlues n(x s ,ti) G Z 0+ . These occupation numbers specify states 

\nUr f.Wnlpsti) 

\n(x„U))= 1 77 JJ , 10)- (Cl) 
y/n{x a ,ti)\ 

From the states of Eq. (Cl), one can form states |n(-,tj)) defined by 

\n(;ti))=U aeZx \n(x a ,ti)) . (C.2) 

Figure C2 shows the input and output arrows for a single node (x s ,ti) of our 
net. This node is assigned a value 

A[n(x s ,t í )|n(·,t í _ 1 )] = (n(x s ,t í )|í]|n(·,í í _i)) , (C.3) 

where Q is a second quantized operator. 

In general, the theory characterized by this QB net is a multi-particle quantum 
field theory. In fact, n(x s ,ti) is the classical field (j)(x s ,U) that is commonly used to 
define Feynman integrals in quantum field theories. In this appendix, we will consider 
a single particle. Thus, for all i, J2 S £Z x n ( x s,ti) = 1 and n(-,ti) = S Xr (-) for some 
r G Z x . For a single particle, the root node (x ,t ) will be taken to have amplitude 
1 for n(x ,t ) = 1 and amplitude for n(x ,t ) = 0. All other nodes will be taken 
to have the value A[n(x s , íi)|n(-, íj-i)] that is specified by the table of Fig.C.3. The 
quantity a s , r in Fig.C.3 is given by 

a s , r = (x s \e^r AtH \x r ) . (CA) 

In this last equation, H is Planck's constant; H is the single particle, first quan- 
tized Hamiltonian; \x r ) for r G Z x are the position eigenvectors, normalized so that 
(x s \x r ) = 5 Sír and J2 r ez x \ x r)(%r\ — 1) where 5 sr is the Kronecker delta. Equation 
(4.1) is trivially satisfied by the table of Fig.C.3. As for Eq.(4.3), it follows from the 
identity 
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\(x 8 \e^-^ H \x )\ 2 = l. (C.5) 



sez x 

2 

For H = f — h V(x, t), one may write 



a S;r = exp 



-íAí 



Ax5(x s — x r ) , (C.6) 



where S(-) is the Dirac delta function. If one were to expand the exponential in the 



last equation, terms in which (-£-) acted only on the delta function would be much 
larger than those in which it acted on V(x s ,t). Thus, one may approximate a S}1 . by 
keeping only those terms in which (-£-) 2 does not act on V(x Sl t). One then gets 



where 



(a s , r ) free e^ AtV ^ , (C.7) 



[a s , r ) free = Ax I 7^e lk{x ° Xr > b ( 2m ) . (C.8) 



oc- 



2tt 



The Gaussian integral in Eq.(C8) is easily performed. One obtains 



where 



a.,r= J— e*^, (C.9) 

V 7T 



A^(f)f(^, (C,0) 



andllTü 



C,r = J (^r) ~V(x s ,t). (C.ll) 

In general, \a s ,r\ « 1- This is why The phase At ^ s ' r of a S:r must not show 
any granularity when x s and x r range over their possible discrete vàlues. Otherwise, 
the outcome of experiments that depended on this phase would depend on Ax. When 
x s goes from x r to x r + Ax, the phase A ^T s,r jumps by A9. If this jump is small, then 



so is |a s ,r| ~- 

Consider Fig.C4, which presents a typical evaluation of the net of Fig.C.l. 
We have blackened those nodes whose state n(x s ,U) is 1. Arrows coming out of 
the blackened nodes are also imagined to lie in state 1. We show only those arrows 
lying in state 1; all other arrows, i.e., those lying in state 0, are not shown. The 
value of each node (the infinitesimal complex amplitude a Str for the blackened nodes 
and approximately 1 for the unblackened ones) has been placed next to it. If one 
connects the blackened nodes one obtains a path. The value of the whole net equals 
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the product of the node vàlues along this path. Thus, the sum over all possible 
evaluations of Fig.C.l with the same final node (at time T) blackened may be viewed 
as a sum over all paths with the same final state. This latter path sum is precisely 
the familiar Feynman path integral that yields the Shroedinger equation for a single 
particle in an arbitrary potential. JÏ5| 

Before concluding this appendix, let us consider the classical limit of the QB 
net of Fig.C.l. 

Recali the method of obtaining the particle's classical equation of motion using 
the action. The action 5'[x(·)] for a path x(t) is 



S[x(-)) = í dtC, (C.12) 



o 



where the Lagrangian £ is 



777 

C = -[x[t)f-V[x[t),t}. (C.13) 

Let 5x(t) be any smooth function of time. For any functional F[x(·)] of a function 
x(t), define 5F by 

5F = F[x(-) + 5x(-)} - F[x(-)} . (CU) 

The equation of motion is obtained by setting SS = to first order in Sx(t) for all 
Sx(-) with Sx(0) = Sx(T) = 0. One gets 

r T ( dV \ 

SS = J dt f -mx(t) - — [x(t),t] J 5x(t) = , (C.15) 

for arbitrary 5x(t), which implies 

dV 

m'x = . (C.16) 

ox 

The classical path x c i(t) is defined as the solution to Eq.(C.16) for given x(0) and 
x{T). The classical action is defined by S c i = S[xd(-)]. When SS = 0, x(t) = x c i(t) 
and 5'[x(·)] = S c i. 

If we restrict x(t) to take on vàlues in the set {x r \r G Z x }, then to the path 
x(-) there corresponds a sequence of integers r(íj) for each i e Z t , where x(ti) = x r ( ti y 
The path x{t) connects a string of nodes (x r ( ti ),ti) (the blackened nodes in Fig.C.4). 
The value of the QB net Fig.C.l is the the product of the vàlues av^),^*^!) of these 
nodes. From Eqs.(C.ll), (C.12) and (C.13), 

%(-)] = E Aí^rfo-i) ■ (C17) 
iez t 

Thus, 
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II a r(ti),r(ti- 

iez t 



-iA9 



7Y 



Nj, 
2 



(C.18) 



Probabilities are calculated by summing the weights Eq.(C.18) for all possible paths 
x(-) with the same x(T) and taking the magnitude squared of the sum. The paths 
which interfere constructively and therefore contribute the most to such probabilities 



are those for which 



<< 1, i.e., those for which the phase change f- accumulated 



over the time interval (0, T) is nearly stationary under deformations 5x(t) of the path. 



<< 1 means SS — > 0. 



In the classical limit, which occurs when % — > 0, 

Note that the classical limit of the QB net Pig.C.l is not its parent CB net. 
With the parent CB net, paths do not interfere, whereas with the QB net in the 
classical limit, there is so much destructive interference between paths that all path 
except the classical one cancel each other out. With the parent CB net, equal weight 
is given to smooth paths like those a classical particle would follow and to jagged 
paths with unbounded variations. With the QB net in the classical limit, paths with 
unbounded variations cancel each other out. 



APPENDIX D. SPIN \ STATES 

Our conventions for spin \ states follows. 

Let ü be a unit 3-dimensional vector characterized by angles (0, (f>) (see Fig.D.l). 

Thus, 



u 



sm V cos ( 
sin 9 sin < 

cos 6 



(D.l) 



Let <7 = (cr x , <Ty, a z ) be the vector of Pauli matrices, where 



1 



-i 



1 



a * = \i o ) > a " = \i o ) ' *' = v o -i 

If and |— u ) are defined by 

<r-u\+ u ) = \+u) , cr-u\- u ) = -\~u) , 

then one can show that 



\+)= CE \ I") 
\+u) [ SE } > I W 




(D.2) 



(D.3) 



(D-4) 



where 
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S = sin^, C = cos^, E = e^. (D.5) 
For example, if 9 = (fi = 0, one gets 

i + ,)=(;), !-•>=(;)■ ( D6 > 

By Eqs.(D.4) and (D.5), if ü and u' are unit vectors characterized by angles 
(9 U , 0u) an d {Qu'i (fiu'), respectively, and if (fi u = <fi u i = 0, then 

(+u'\+u) = (-u'\-u) = cos ( ^— — ) , (D.7) 



(+«'|-«) = -K'I+u) = «in y 2 J • (D-8) 
For any unit vector ü, the antisymmetric state for two partides 1 and 2 is 

IC nt ) = -^[|+ u )i|- u ) 2 -|- u )i|+«) 2 ] • (D.9) 

(IVC 1 *) nas zero total angular momentum and is thus also called the singlet state). 
By expressing the states |±„) on the right hand side of Eq.(D.9) in terms of |± 2 ), it 
is easy to show that \ifi'^ lt ) is invariant under rotations, i.e., it is independent of ü. 

So far, our formulation has been a first quantized one. In a second quantized 
formulation, one defines annihilation operators a ua , for any unit vector ü and for 
a G {+, — } • These annihilation operators must satisfy 

[cW,a u / CT /]+ = , (D. 10) 

where, for any two operators A and B, [A, B] + = AB + BA. States in the first 
quantized formulation are mapped in a 1-1 fashion into states in the second quantized 
formulation. For example, 

|+ u > - aí + \0) , (D.12) 

|-«>-a£-|0), (D.13) 

\r nt ) ^ aí + al-\0) ■ (D.14) 
From Eqs.(D.4), (D.12) and (D.13), it follows that 
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±)-(-sé- S c E E )(ÍV 



"z— 



(To check this last equation, just apply |0) from the right to both sides of the equa- 
tion.) 
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Figure CAPTIONS: 

Fig.l Node labelled by the random variable Xj. 
Fig.2 Possible CB nets with 2 nodes x and y. 
Fig.3 Possible CB nets with 3 nodes x, y and z. 
Fig.4 Three-node cycle. 

Fig.5 Proof that the joint probability represented by Fig.3d adds up to one (i.e., 
satisfies Eq.(2.4) ). Summation over the states of an arrow is indicated by 
giving the arrow a double shaft. 

Fig.6 Figure 6a shows the fully connected four node graph with its chronological 
labelling. By deforming Fig.6a into a topologically equivalent diagram, one 
obtains a more "stylized" version of the same thing, Fig.6b. 

Fig.7 Chapman-Kolgomorov equation for a Markov chain (x,y,z). 

Fig.8 Generalization of Fig.7 to arbitrary random variables x,y,z that don't neces- 
sarily form a Markov chain. 

Fig.9 CB net version of an AND gate. 

Fig.10 z is a SUM node. 

Fig.ll z is an IF-THEN node. 

Fig.12 CB net for the Clauser-Horne experiment. This figure really represents 4 nets, 
oneforeachofthefollowingpossibilities: (0i,0 2 ) = (A, B), (A, B r ), (A! , B), (A\ B'). 

Fig.13 CB net for an experiment like the usual Clauser-Horne experiment, except 
that here the measurement types Q\ and 82 are selected at random. 

Fig.14 CB nets associated with a random walk. Net Fig.l4c is coarser than net 
Fig.l4b, which is coarser than net Fig.l4a. 

Fig.15 Root node representing initial wavefunction. 

Fig.l6a Marginalizer node. 

Fig.l6b Phase Shifter node. 

Fig.17 Nodes for Stern-Gerlach magnet with (a) one, and (b) two incoming modes. 

Fig.18 QB net for an experiment with 2 Stern-Gerlach magnets. (triangle=wavefunction, 
black-filled circle=marginalizer, white-filled circle=magnet). 
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19 QB net for an experiment with 2 Stern-Gerlach magnets. (triangle=wavefunction, 
black-filled circle=marginalizer, white-filled circle=magnet). 

20 Various sets of evidence that were considered for the nets Figs.18 and 19 with 
two Stern-Gerlach magnets. 

21 One-hypothesis probabilities, for the net of Fig. 19, for evidence-cases 1, 2, 4, 
10 and 12. Columns A to D refer to the parent CB net and columns F to I to 
the QB net. 

22 Two-hypothesis probabilities, for the net of Fig. 19, for evidence-cases 1, 2, 4, 
10 and 12. Columns A to F refer to the parent CB net and columns H to M to 
the QB net. 



23 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

24 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

25 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

26 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

27 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

28 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 

29 QB net for an experiment with 3 Stern-Gerlach magnets. 
black-filled circle=marginalizer, white-filled circle=magnet 



triangle=wavefunction, 
triangle=wavefunction, 
triangle=wavefunction, 
triangle=wavefunction, 
triangle=wavefunction, 
triangle=wavefunction, 
triangle=wavefunction, 



30 Various sets of evidence that were considered for the nets Figs.23 to 29 with 
three Stern-Gerlach magnets. 

A.l Graphical proof of an identity about the intersection of two direct product 
sets. 

Cl "The fabric of spacetime". QB net that yields the Feynman path integral 
for a single mass particle in an arbitrary potential. 

C.2 The input and output arrows for a single node (x s ,ti) of the net Fig.C.l. 
C.3 In the net Fig.C.l, the value A[n(x s: ti)\n{-, íj-i)] of the node (x s ,ti). 
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Fig.C.4 A typical evaluation of the net of Fig.C.l. We have blackened those nodes 
whose state n(x s ,ti) is 1. Arrows coming out of the blackened nodes are also 
imagined to lie in state 1. We show only those arrows lying in state 1; all other 
arrows, i.e., those lying in state 0, are not shown. The value of each node (the 
infinitesimal complex amplitude a s . r for the blackened nodes and approximately 
1 for the unblackened ones) has been placed next to it. 

Fig.D.l Definition of the angles 9 and <fi used in specifying the unit vector ü. 
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